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Abstract 

Having first performed a Monte Carlo simulation to justify the analysis technique 
to be used, we then analyze the Bose-Einstein correlation functions from CERN exper- 
iment NA44 in the context of the core-halo model. Although experimental resolution 
and error bar distribution prevents a direct observation of the halo structure, the 
values for the core radius and the core fraction of pions can be obtained in a straight- 
forward manner. These are found to be independent of the structure of the correlation 
function at small relative momenta of Q < 50 MeV. Hence, the cj meson decay prod- 
ucts do not distort the Gaussian shape of the Bose-Einstein correlation function in 
S + Pb^TT + TT + X reaction at CERN SPS. As we find that the "model-independent" 
HBT radii yield results that are quantitatively as well as qualitatively unreliable for 
systems with long-lived resonances, we present their corrected form that applies for 
correlation functions with A(K) < 1. 







1 Introduction 



High energy heavy ion colhsions are providing a reproducible experimental environment to 
study physics at extremely high energy densities in relatively large volumes. A tool to access 
the space-time characteristics of such systems is the technique known as intensity interferom- 
etry. If the observed particles are bosons, intensity interferometry is frequently referred to as 
HBT effect or Bose-Einstein correlations ^, ^ . These intensity correlations appear due 
to the enhanced (decreased) likelihood that a boson (fermion) is produced in a quantal state 
that is close in phase-space to a state already occupied by another boson (fermion). Intensity 
interferometry is a fastly developing tool for studying spatial and temporal extension of hot 
and dense strongly interacting systems on the 10~^^ m and 10"^'^ s scales. 

Accumulating evidence indicates, that high energy heavy ion reactions, performed cur- 
rently at CERN SPS and Brookhaven AGS, create composite sources of particles, that can 
be approximately divided into two parts: a central core, where the time evolution of the 
matter can be thought as a violent cascading of binary collisions or hydrodynamical behav- 
ior of strongly interacting matter, and a halo of long-lived resonances, that are created in 
the hot and dense medium but decay far outside the core. In case of pions, the resonance 
halo undoubtedly contains the decays of where r = 8.9 ■ 10~^ s, the decays of the t] and 
rj' resonances that have a very narrow decay widths of = 1.20 keV and F^/ = 0.2MeV. 
The lifetimes of these resonances correspond to r = h/T = 164417 fm/c and 986.5 fm/c, 
respectively. The next long-lived resonance which is expected to be produced in significant 
numbers in high energy collisions is the u meson with a full width of F^ = 8.4 MeV and 
a corresponding life-time of r^^ = 23.4 fm/c. In the literature, one may observe a scientific 
debate weather the decay products of the a; mesons contribute to the halo or not. 

After a brief summary of the statements made about the contribution of the omega 
resonance to the BECF-s, we attempt to conclude this debate in case of the S+Pb tt+tt+X 
reaction measured by the NA44 Collaboration at CERN SPS 0, by presenting an analysis 
of this data set along the lines suggested in ref. [§. In case of other reactions or other 
experimental resolutions/techniques, we emphasize that the role of the u meson should be 
investigated again, in a manner similar to the analysis presented in the forthcoming sections. 

The main purpose of the present paper is to test, whether the core/halo model inter- 
pretation of the two-particle correlation function is indeed applicable to the above cited 
NA44 data, or not, with the help of a novel method suggested in ref. P|. If the core/halo 
interpretation of the measured correlation function is meaningful, the parameterization of 
the correlation data should be independent of the data points at the lowest values of the 
relative momentum, hence these data points could be deleted without significantly changing 
the fitted parameters of the correlation functions. We check the validity of this hypothesis 
by gradually removing the data points at smallest relative momentum from the NA44 S + Pb 
data set. To our best knowledge, this is the first time that such kind of analysis is performed 
on measured correlation data, in order to verify the applicability of the core/halo model. 
Hence, we must establish the limits of the technique to be used before attempting the data 
analysis. 

The paper is structured as follows: In Section ^ we summarize the basics of Bose-Einstein 
correlations and the assumptions that are introduced to obtain the core/halo model. We also 
summarize here those various theoretical considerations and conclusions about the influence 
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of the uj meson on the two-pion correlation function, that were achieved in earher numerical 
investigations, reviewing the status of this scientific debate. 

In Section |^, we describe a Monte-Carlo procedure that was used to estimate the limits 
of the procedure that we utilized to verify the applicability of the core/halo picture to NA44 
data in the subsequent Section ^. Finally, we summarize and conclude. Conceptually, our 
study is also related to a recently discussed "model-independent" Gaussian parameterization 
of the correlation functions. In Appendix A we evaluate the Gaussian model-independent 
radii in the core/halo model and compare the resulting "model- independent" correlation 
functions to the prediction of the core-halo model. Then we give a re-formulation of the 
"model-independent" parameterization of the correlation functions that can be applied to 
core/halo type of models as well as systems without halo and clarify the model assumptions 
that have to be made in order to achieve this "model-independent" result. 

In the foUowings, we use natural units, ^ = c = 1, to simplify the notation. 



2 Correlation Functions and the cj-Puzzle 

Let us briefly recapitulate the theoretical description of the two-particle Bose-Einstein cor- 
relation functions following the lines of refs. 0, |^ . 

Particle emission is characterized by the single-particle Wigner- function S{x,p). Here 
X = {t,r) denotes a point in space-time and p = [E, p ) a point in momentum-space. The 
emitted particles are on mass-shell, m = \/E'^ — \)'^ . 

An auxiliary quantity can be introduced as 



^(AA;, K) = J (fx S{x, K) exp{iAk ■ x) 



where Ak = pi — P2, K = {pi + P2)/2 and Ak ■ x denotes an inner product. 
The invariant momentum distribution can be expressed as 

E— = N,ip) = SiAk = 0,K = p). (2) 
dp 

In the present paper, we utilize the so called hydrodynamical normalization of the Wigner- 
functions. 



^d'xS{x,p) = J -fN,{p) = (n), (3) 

where (n) is the mean multiplicity. 

The two-particle correlation function can be written as 

g.(P.,P.) = ./^'P';/^' > = 1+ y:'-:^:' ^ ^ i+ . (4) 
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e.g. in ref. [§, |^. Interestingly, this formalism can be introduced not only when multi-particle 
symmetrizations are negligible, but even when they are fully taken into account |jlO|- For 
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the most recent derivation, where full multi-particle symmetrization effects were included, 
see ref. [|10 . 

The last approximation in eq. (Q) can be estimated to give correct results within 5 % 
error [O]. We have neglected here final state interactions, and a completely chaotic particle 



emission is assumed. 

The second and the higher order Bose-Einstein correlation functions. 



Cn(Pl,P2, •••,P^ 



are given 



iVi(pi)iVi(p2)...iVi(p„) 
in terms of the Fourier-transformed Wigner-f unctions as 



(5) 
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where cr"' stands for the set of permutations of n indexes, 
S{i,ai) = S'(Ki,^^, Aki,^J, and s{i,ai 

Pi + Pa, 



and Ak,- , 



S{i,cri) 
S{i,i) 

Pi - Pa,- 



(7) 
(8) 



Note the difference between a", which stands for a set of permutations and Uj (subscript i), 
which stands for the permuted value of index z in a given permutation from the set a". 
Let us recapitulate the Assumptions of the core/halo model from ref. 0: 
Assumption 0: The emission function does not have a no-scale, power-law like structure. 



This possibility was discussed and related to intermittency in ref. |13 



Assumption 1: The bosons are emitted either from a centra/ part or from the surrounding 
halo. Their emission functions are indicated by Sc{x,p) and Sh{x,p), respectively. According 
to this assumption, the complete emission function can be written as 



S{x,p) = Sc{x,p) + Sh{x,p), 



(9) 



using the hydrodynamic normalization of the Wigner functions. 

Assumption 2: We assume that the emission function which characterizes the halo changes 
on a scale Rh which is larger than Rmax ~ ^IQmini the maximum length-scale resolvable 
by the intensity interferometry microscope. However, the smaller central part of size Rc is 
assumed to be resolvable, Rn > Rmax > Rc- This inequality is assumed to be satisfied by 
all characteristic scales in the halo and in the central part, e.g. in case the side, out or 
longitudinal components [|l^, of the correlation function are not identical. 

Assumption 3: The momentum-dependent core fraction /c(^) = Nc{pi)/Ni{^i) varies 
slowly on the relative momentum scale given by the correlator of the core Sc(l; 2)sc(2, 1). 

Let us also recapitulate the normalization conditions: 



j d^x-^Sc{x,p) = (n)c, and J d'^x-^Sh{x,p) = {n)h 



dp 



(10) 
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where the subscripts c, h refer to the contribution from the central core and from the halo, 
respectively. 

Note that ref also utilized the above Assumptions, however, in a modified form, 
that was based on Wigner-functions normalized to 1, and a core-fraction fc = {n)c/ (n) 
was introduced, while in the present paper a momentum- dependent core fraction /c(p) = 
Nc{p)/N[p) is utilized. One finds that 

N,{p) = N,{p) + NM, and {n) = {n),+ {n)f,. (11) 
According to ref. 0, the general expression for C„(l, n) reads in the core/halo model as 

C„(l,...,n) = 1 + E E Efl/c(^'^)5cfevJ- (12) 

j=2h,...,ij=l p" k=l 

where J2' refers to a summation over different values of the indices (i.e. ij = ik type of 
terms are excluded). For the two-particle correlation function, the above equation takes a 
particularly simple form: 

C2(l,2) = l + /,(l)/,(2)§e(l,2)5e(2,l). (13) 

With the help of Assumption 3, the core/halo model thus predicts the following form for the 
two-particle correlation function: 

|2 



r./Ai ^ 1 , \ ^ I 5'c(Aki2,Ki2) I i , \ ^ | 5'c(Aki2, K12) 

C(Aki2,Ki2) = 1 A4Ki2)^7;^ ^^7;^ ^, ^l+A*(Ki2) , ~ 



^,(0,Pi)^,(0,P2) |^c(0,Ki2 

where one introduces an effective intercept parameter as 



(14) 



A.(K) = [iV,(K)/A^l(K)]^ (15) 

As emphasized in Ref. 0, this effective intercept parameter ( which is not the same as the 
exact intercept parameter, A^,. = 1 at Q = MeV) shall in general depend on the mean 
momentum of the observed boson pair, which within the errors of Qmin coincides with any 
of the on-shell four-momentum pi or p2. 

Thus one obtains the core/halo interpretation of the two-particle correlation function: 
The measured part of the BECF picks up an effective, momentum dependent intercept 
parameter A*(K), that depends on the mean momentum only and which can be used as 
a tool to measure the momentum dependence of the fraction of particles that are emitted 
from the core. On the other hand, the relative momentum dependence of the BECF-s 
carries the information on the core in this picture. Note, that in high energy heavy ion 
collisions the momentum dependence of the A*(K) parameter is very weak, actually, within 
the errors A*(K) is constant for the NA44 data analyzed in ref. 0. However, the validity 
of Assumption 3 has to be checked experimentally for each data set, by determining the 
momentum dependence of the A(K) parameter of the two-particle correlation function. 

At this point, we emphasize that non-Gaussian correlation functions with < A=k < 1 are 
very well possible within the core/halo picture, as discussed e.g. in refs. [^, [1^. We shall 
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discuss in Appendix A under what conditions can the core correlator have an approximate 
Gaussian shape. 

Since the core-halo model neglects possible partial coherence, motivated by the success of 
fully chaotic Monte-Carlo simulations of high energy heavy ion collisions like RQMD |jT^, the 



exact intercept parameter is = 1 and the value of the exact intercept of the two-particle 
correlation function is always 1 + A^; = 2 in the core/halo picture. The effect of long-lived 
resonances is to create an unresolvable, narrow peak in the region < Q < Qmin from the 
interference of the particles of (/i, h) and c) type, i.e. from those pairs which have at least 
one member of the pair from the core. This is illustrated in Figure |l|, taken from Ref. 

Another very important general observation is that the region with small Q has a different 
characteristic structure from that oi Q > Qmin region. Thus any analytical expansions of 
the Bose-Einstein correlation functions around the (5 = point become both qualitatively 
and quantitatively unreliable for the core-halo type of systems, characterized by < A* < 1. 
This point is further elaborated in the Appendix, let us return now to the detailed discussion 
of the core-halo model using an example. 

Note, that in principle the core as well as the halo part of the emission function could 
be decomposed into more detailed contributions. We shall argue below, that in case of 
5* + Ph reactions at 200 AGeV at CERN SPS, in the NA44 acceptance, one can separate the 
contribution of various long-lived resonances as 

S,{x,p) = Sl[\x,p) and Ar,(p) = ^ Nl'\p). (16) 

r=uj,r},r]' ,Kg r=uj,r},r]' ,Kg 

For a general consideration, e.g. the ones discussed in ref. 0, the question whether 
the u decay products contribute to the halo or not, is essentially indifferent. However, 
when a data analysis is performed, this question becomes important both qualitatively and 
quantitatively, influence of the u meson on the shape and the intercept parameters of the 



two-particle correlation functions is suggested in various theoretical papers, see e.g. refs. ||T5 
|T3, 0, 0, |2T|, ^, ^ 



Let us briefly summarize the various statements made in this scientific debate. Essentially, 
three different type of statements were made regarding the influence of decay products of 
the u meson on the short range part of the two-particle correlation function. Although we 
cannot repeat the fine-prints and all the reservations that were made in earlier papers on this 
point, we think that a rough (and probably incomplete) summary can be made as follows: 

A The effect of the u decay products is clearly visible on the two-particle BECF. These decay 
products are mainly responsible for the deviation of the BECF from a Gaussian shape 
and they influence the intercept parameter also in a non-trivial manner. Essentially, 
the u decays contribute to the resolvable, core part of the source. 

B The decay products of the u meson should be considered as an intermediate case. Their 
contribution should be partially visible, resulting in a deviation from the Gaussian 
shape of the BECF. 

C At the current level of experimental techniques, Qmin ~ 5 — 10 MeV and B{Q)max > 25 
MeV, the uj decay products can be taken as part of the halo. Their contribution to the 
two-particle BECF cannot he resolved experimentally in case of NA44 data. 
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In the core/halo picture, the points A - C correspond to the following conclusions: In 
case of A, the uj decay products belong to the core. In case of B, the core/halo model is 
not applicable to simplify the theoretical treatment of the two-particle BECF-s. This case 
corresponds to a really complicated situation. Finally, in case of C, the u decay products 
contribute to the halo and the theoretical description of the correlation function can be 
simplified substantially. 

Let us also give a (probably incomplete) list of papers where these statements were made. 
Conclusion A has been reached in ref. |T^, in a Monte-Carlo simulation with the help of 
SPACER, and later in ref. |^ assuming 5 MeV Qmin resolution for Ph + Ph reactions at 



CERN SPS in a 3d hydrodynamical simulation with the code HYLANDER. Conclusion B 
was reached in refs. pO|, ^ , using phenomenological parameterization of a hydrody- 
namically evolving core with thermally populated resonance production. In all these studies, 
full resonance decay kinematics was included and resonance production was described either 
in a thermal manner or in a non-equilibrium process with the help of re-scattering. All of 
these studies share in common that they evaluated the correlation function as a mathemati- 
cal function but they did not consider in detail the error bar distribution on the theoretical 
curve. 

Finally, in refs. @, ^, ^ 26 1, conclusion C was reached based on considerations related 



to the NA44 detector resolution and the surprisingly Gaussian shape of the measured NA44 
correlation functions. 

Another interesting study was performed recently by Padula and Gyulassy [0, that 
included a simulation of binning and detector resolution effects for AGS energies with the 
help of the inside-out side cascade code CERES, with proper resonance decay kinematics, to 
conclude that 

D Preliminary data seem to rule out dynamical models with significant uj and rj resonance 
fraction yields in case of the E802 measurement of Si + Au data at 14.6 AGeV at 
Brookhaven AGS. 

We do restrict our study to CERN SPS reactions, and published NA44 data. In the 
next section, we perform a Monte-Carlo simulation to make a connection between source 
parameters and the parameters of the fitted correlation functions before attempting the 
data analysis. 



3 Monte Carlo Simulations 

For the sake of simplicity, one may simplify the problem by fixing the mean momentum of 
the observed particle pairs and introduce the source density for those particles which are 
emitted with the given mean momentum as p-k{x) = S{x,K.). The index K will thus be 
suppressed in the forthcoming, but implicitly we shall assume that the analysis is performed 
at a fixed value of K. Also, in the applications, we shall analyze NA44 data where the mean 
momentum of the pairs is restricted to a certain range. We have had access only to the one 
dimensional slices of a three-dimensional NA44 data set, along the main axis in the LCMS 
frame |]T5|. These data were taken at a fixed mean momentum. In effect, we thus analyzed 
three different one dimensional data set instead of analyzing a single three dimensional one. 
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As the quality of the pubhshed data is expected to be improved in the near future, we hope 
that this hmitation will disappear and in a future analysis one may study the 3d correlation 
functions at a fixed value of K directly in a similar manner as presented below. 

The limitation that we had access to Id slices of data yields an advantage as well, since it 
is enough for the present paper to formulate the procedure for one dimensional distributions 
only. The generalization of the method for multi-dimensional distributions is trivial hence 
it will be omitted. 

In the core/halo picture, the source density is written as 

p{x) = fcPcix) + (1 - fc)Ph{x) (17) 

where fc from now on represents the fraction of pions produced in the core at the given 
mean momentum K. In order to make predictions for the model, forms for Pc{x) and Ph{x) 
must be assumed. With Gaussian assumptions for both, pi{x) = ^^^^ e"^^/^'^? , the resulting 
correlation function is 

C2{Q) = 1 + /.^e-^?«^ + (1 - f.re-^l^' + 2U1 - /Je-^(^?+^^)«^ (18) 

Examination of this expression reveals Gaussian terms for the core, the halo, and an 
'interference' term. Figure ^ shows the relative size of these terms, as a function of some 
momentum variable Q, for typical values of fc, Rc, and Rh- The effect of the halo is to 
introduce a sharp peak in the correlation function at low values of Q. This narrow peak 
in Q-space is indicative of a large length-scale in x-space. Note that the model predicts 
C2{Q = 0) = 2 - there is no need to assume any coherent boson production which would 
reduce the intercept. Figure ^ illustrates the effect of the fc parameter on the correlation 
function. 

Since we are more interested in the radius parameters of the core than the halo, the 
easiest way to extract the proper core parameters is to remove the low Q data points so that 
the effect of the halo is negligible, a procedure suggested e.g. in ref . . This simple method 
minimizes the effect of the Gaussian assumption for the shape of the halo structure - as 
long as the length-scale of the halo is large compared to h/Qmin, removing the low Q points 
should remove the halo's effects. Fitting the remaining data with a standard Gaussian form 
should produce R = Rc and A = 

However, given the finite experimental resolution of ~10 MeV/c, the halo may not be 
visible at all. If this is the case, the extracted A and R parameters still have an interpretation 
as the fraction of core pions (squared) and the core radius parameter, respectively. If they do 
not change as Qmin, the minimum value of Q included in the data, is increased, it will indicate 
that the correlation function really is Gaussian. This would contradict the prediction that 
some resonances distort the Gaussian shape of the correlation function. 

In either case, we must establish the limits of the procedure before attempting the analysis 
of real data. There must be some maximum number of data points that can be removed from 
a real correlation function while still being able to extract valid parameters. If that number 
is too small, then the idea of removing low Q points in order to study the correlation function 
of the core is invalid - if the length-scales of the core and the halo cannot be separated, then 
the fitted parameters of the truncated correlation function become dependent on Qmin, the 
size of the excluded region. 
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The core-halo model takes advantage of finite experimental resolution - as resolution 
increases, and more of the low Q structure of C2{Q) is revealed, perhaps a more detailed 
model will be necessary. But for now, this reasonable model imitates the effects of resonance 
production and can attribute physical significance to the R and A parameters extracted from 
a correlation function. 



3.1 The acceptance-rejection method 

What is required is a method of producing correlation functions where the underlying pa- 
rameters are known in advance, so they can be compared to the extracted parameters. This 
will be done by simulating the actual distribution A{Q) and the background distribution 
B{Q) that a real experiment uses to produce C2{Q) = A{Q) / B{Q). 

It is necessary to assume a form for A{Q) and B{Q) in order to do this. The background 
distribution B{Q) is different for every experiment, so one must parameterize it appropriately 
each time. The A{Q) distribution will not be parameterized directly; rather, we will assume 
a form for C2{Q), and then simulate the distribution A{Q) = B{Q) ■ C2{Q)- In this way, 
we can have control over the parameters that are going into the simulation. The finite 
number of iterations performed in the simulation ensure that these parent distributions are 
never reproduced exactly, so the assumed form for C2{Q) is not returned directly. Binning 
effects are automatically included and, most importantly, the error-bar distribution on the 
correlation function will follow the experimental one at least when statistical errors are 
considered. Fits shall be done with optimizing /NDF hence the error bar distribution 
plays an important role when determining the best fitted values and their errors. 

We would like to emphasize that the simulation of the A{Q) and B{Q) distributions is 
essential when discussing possible resonance effects on the two-particle correlation function, 
since the resonance decay products are expected to influence the shape of the Bose-Einstein 
correlation function in the Q < 50 MeV region |]15|, |T9|, ^, ^ |23| and they also contribute 



to the fltted value of the radius and the intercept parameters. However, the fitted values 
are strongly infiuenced by the distribution of the errors on the data points, since the fit 
selects to reproduce the best that part of the correlation function, which has the smallest 
errors. Unfortunately, very few experiments decided to publish their actual and background 
distributions, but some |28|, ^ |3y, Q did this important step. Even fewer theoretical 



models tried to reproduce the distribution of the statistical errors on the data as arising 
from the number of actual and background pairs as a function of the relative momentum. 
In fact, none of the theoretical simulations that reached a conclusion belonging to type A 
or B included this important step. 

Unfortunately, it is not possible to sample from the complicated distributions A{Q) and 
B{Q) directly. Instead, we will use the acceptance-rejection method, which works as follows: 
Let the known distribution to be sampled be denoted by p{Q), that will be chosen as A{Q) 
or B{Q) in the subsequent parts. The probability distribution p{Q) should be normalized, 
/ p{Q) dQ = 1, over the range of interest. Find a comparison function, f{Q), which has the 
following properties: 

1. Let f{Q) satisfy f{Q) > p{Q) for all Q. 

2. There exists a closed form for a random Q point sampled from f{Q). 
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In the interest of computational efficiency, the comparison function f{Q) should be chosen 
to be as close to p{Q) as possible, while still satisfying f{Q) > p{Q), in order to minimize 
the number of rejected points. A random Q point, Qi, is sampled from the comparison 
function f{Q). Then the value of the parent distribution p{Qi) is calculated. The ratio 
R = p{Qi)/ f{Qi), which is always less than 1, is then compared to another uniform random 
number U'. If U' > R, the point Xi is rejected and the process starts over again. If U' < R, 
Qi is accepted. In this way, the distribution p{Q) is sampled. 



3.2 Simulating a one dimensional correlation function 

If we assume the halo form for C2{Q), eq. the two distributions to be simulated are B{Q) 
and A{Q). For the CERN experiment NA44, studying S+Pb collisions, the shape of B{Q) 
is given approximately by ^2 

B{Q) = g3e-3.6Q-^ (19) 
hence the A{Q) distribution is sampled as 

AiQ) = 5(g)(l + /^-^^«^ + (1 - f.re-^l'^' + 2/,(l - /,)e-^(«^+^^^)«^). (20) 

For these one-dimensional simulations, Q can be any relative momentum variable - Qmv, 
Qside, Qout, etc. The above expression for B{Q) was given for Qinv, but it is valid for 
any one dimension of Q as long as the other components are small. The parameters for a 
Lorentzian comparison function are adjusted by hand until the Lorentzian satisfies condition 
1) reasonably well. The Lorentzian comparison function was sampled as suggested in ref. |P3|| . 
With the help of Lorentzian comparison functions, the distributions A{Q) and B{Q) were 
sampled as described above. As each point was sampled from either A{Q) or B{Q), it is 
binned. The bin size of 10 MeV/c is chosen to reflect experimental situation for the multi- 
dimensional NA44 data analysis. The value of Q for each bin is calculated by averaging 
all entries, from both A{Q) and B{Q), that fall in that bin. Figure ^ shows the results 
of the simulations of A{Q), B{Q), and the resulting C2{Q)- The number of iterations for 
the simulation is chosen to give reasonable error bars on the final C2{Q)- The number of 
iterations must be higher for A{Q), however, since A{Q) > B{Q) for all Q. The areas under 
the curves are used to calculate the proper ratio of iterations. 

The error bars are calculated as follows. For A{Q) and B{Q), we assume that the number 
of entries in each bin is Poisson distributed, so that the error for each bin is just the square 
root of the number of bin entries. The errors for C2{Q) are then just combined in the usual 
manner. 

Now we can fit the simulated correlation function with a standard Gaussian shape, 
C2{Q) = 1 + Xe~^^'^^ , reflecting experimental analysis techniques, and observe the behavior 
of the parameters as Qmin is increased. Figure ^ shows the value of the i?* and A=k parameters 
as low Q data points are removed. The solid lines represent the input values, Rc = 4.0 fm 
and = 0.5625. Examining these results indicates that reliable values for the parameters 
can be extracted, i.e. R^, = Rc within errors and A.,, = within errors, when Qmin is varied 
in the range of 10 to 50 MeV/c. We can also read off from Figure |^ that the fitted radius 
and intercept values are essentially insensitive to the exact value of Qmin if this is varied in 
the above range. 
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4 Core/halo model analysis of NA44 data 



We now apply the data chopping analysis, tested in the earlier section, to real experimental 
data. CERN experiment NA44 has been running since 1991, and has an extensive collection 



of correlation function data [29, 30, 34, 35, 31, pi. We analyze their S+Pb data. The data 



set is three-dimensional - that is, it covers all of the (out, side, long) space, but we only 
have access to slices of the data along the three axes. So in effect we are analyzing three 
one-dimensional slices of the three-dimensional Bose-Einstein correlation functions, for each 
published data set. Values of the two-particle correlation functions and their errors were 
read directly off the published figures by eye, as accurately as possible. NA44 has analyzed 
S+Pb data for both 7r+ and . The 7r+ data were taken with two different magnetic 
settings - 'high pt data is characterized by < pt >= 450MeV^/c, and the 'low pt data has 
< Pt >= 150MeV/c. This was done in an attempt to observe some detailed dynamics of the 
system. 

Thus we have analyzed the slices of the three-dimensional NA44 Bose-Einstein correlation 
function for the \ow-pt and high-p^ pion sample and for the kaon sample in S" + Pb reactions 
at CERN SPS, data were from refs. 0. Unfortunately, we had no access to the three- 
dimensional distributions either for NA44 measurements or for other data. Ideally, the 
analysis presented below should have been performed on the (unpublished) three-dimensional 
distributions. However, our final result indicates that the extracted values for the intercept 
parameters are within errors similar for each projection of the NA44 sample, which justifies 
the utilization of the projections in this particular case. 

Our analysis results are shown in the Figures |^-|TT|. Note that the left-most data point 
in each figure represents the value of that particular parameter when none of the data 
points are removed. Also, with three separate slices along the three main axis instead of one 
three-dimensional set, one can extract three values for the A parameter. The NA44 analysis 
produces just one A, since they fit to 

C2{Q) = 1 + Ae~'^°"*^°"'~*^^*'*'=^^"^<=~'^^°"9'^^°"9 (21) 

The parameters show no significant, systematic change as Qmin is increased. This con- 
firms that the shape of the correlation function really is Gaussian. In the context of the 
core-halo model, this suggests that long-lived resonances are not resolved by NA44, neither 
in case of pions, nor for kaons, see Figures (|6|-p!l|). 

In order to extract values for R and fc from this analysis, we fit the previous data with 
a constant (see Table |l]). The data up to Qmin = 40 MeV/c are chosen, since this has been 
shown to be a reliable range. 

These results show several things. As expected, the kaon core fraction is larger than 
either pion fraction. Further, the trend in the core radii is Rkaon < Rhighpt < Riowpt- The 
interpretation of this result is not a trivial task One would first imagine that the 



decrease of these radii follows the trend of increasing resonance influence, suggesting that 
short-lived resonances can also affect core radii. It is also possible that kaons decouple from 
the interaction region sooner than pions do. Secondly, the core fractions are similar for 
the high and low pt data, since the intercept parameter is within errors independent of the 
transversal mass of the particles. One might argue that the constancy of the A parameter 
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parameter 


high Pt 7f^ 


low Pt 71^ 


K+ 


Rc out (fm) 

Rc side (fm) 
Rc long (fm) 
fc out 
fc side 
fc long 


2.92 ± .13 
2.90 ± .18 
3.31 ± .16 
.704 ± .012 
.735 ± .021 
.738 ± .014 


4.29 ± .13 
4.24 ± .26 
5.43 ± .30 
.725 ± .011 
.647 ± .021 
.724 ± .015 


2.54 ± .18 
2.22 ± .19 
2.67 ± .22 
.802 ± .027 
.736 ± .033 
.789 ± .029 



Table 1: Extracted values for Rc and fc from NA44 S+Pb data. 

implies that any increase in resonance contribution in the low pt data takes place in the 
short-lived resonances. 



Name 


m (MeV) 


r 


Decay 


fpritiof 


fRQMD 


direct vr 


140 


7.804 m 




0.19 


0.33 


P 


770 


1.3 fm 


TTTT 


0.40 


0.26 


A 


1232 


1.64 fm 


Ntt 


0.06 


0.12 


K* 


892 


3.94 fm 


Ktt 


0.09 


0.07 


S(1385) 


1385 


5.5 fm 




0.01 


0.02 


UJ 


783 


23.4 fm 


TTTTTT 


0.16 


0.07 


Tjl 


958 


982 fm 


rjmT 


0.02 


0.02 


7] 


548 


1.64 X 10^ fm 


TTTTTT 


0.04 


0.03 




498 


2.7 cm 


TTTT 


0.03 


0.07 



Table 2: Sources of pions in relativistic heavy ion collisions, taken from Ref. |2^ 



This interpretation, however, is not likely since hydrodynamical expansion may also result 



in a strong transversal mass dependence of the measured radius parameters |^, even an 
approximate rrit scaling can be obtained in this manner in a certain limiting case. Further, a 
comparison of Tables |I] and ^ indicates that the measured core fractions are in the vicinity 
of the core fractions of Fritiof and RQMD if u, rj and rj' are taken as unresolved long lived 
resonances, as indicated by Table ^ 



fc,NAu (high Pt TT+ ) 


fc,NAU (low Pt TT+ ) 


fc, Fritiof 


fc,RQMD 


0.71 ± 0.01 


0.71 ± 0.01 


0.75 


0.80 



Table 3: Core fractions fc from NA44 S+Pb data compared to core fractions from RQMD 
and Fritiof, when u, t] and t]' are taken as unresolved long-lived resonances. The core 
fractions for the NA44 low pt and high pt sample are obtained from a simultaneous fit to 
{fc,out, fc,side, fc,iong) by a coustaut value, and the result is rounded to two decimal digits. 

Since both uj and r] decay predominantly into three pions, in a narrow region of phase 
space characterized by m — 3m7r, these resonance decays should produce predominantly 
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low-pt pions. Similar conclusion holds also for the t]' resonance where the phase space is 
characterized by m^/ — 2m7r — m^, which is also rather small. Thus a resonance interpretation 
of the rrit dependence of the radius parameters seems to contradict to the observed constancy 
of the A parameters. Collective effects, like three dimensional hydrodynamic expansion 
and rrit dependent volume factors that enhance the direct production of pions at low pt as 
compared to the direct production of heavier resonances were discussed as a possible solution 
of this puzzle in Refs. [|16|, |, H. 

In the future, we hope that this type of analysis will be applied to full three-dimensional 
data sets by various experimental groups. This would be more significant for several reasons. 
First, several recent papers |3^ have suggested the presence of cross-terms in the 3-D 
correlation function, something which cannot be detected with just slices of data along the 
three axes. Secondly, the 1-D slices we have used may cover up some of the \ow-Q detail, 
since they integrate over small regions of the other two dimensions (typically < 40 MeV/c). 

Looking farther ahead, the core-halo model could be useful in the next generation of 
heavy ion experiments. Physicists at RHIC or at CERN experiment NA49 expect to have a 
Q resolution of < 5 MeV/c, which could be sufficient to see the halo directly. 

5 Conclusions 

In summary, an outline of Bose-Einstein correlation analysis for heavy ion collisions has been 
presented. In contrast to theoretical predictions of a complicated correlation function due to 
the effect of resonance decays, experimental correlation functions seem to be characterized 
by just a single intercept parameter A and a set of Gaussian radius parameters, Rside, Rout, 
Riong (and Rout-long)- As an attempt to incorporate resonance effects into a simple correlation 
model, core and resonance halo sources are dealt with separately, each characterized by a 
Gaussian width parameter along with a core fraction parameter. Data from NA44 do indicate 
indirectly the presence of such a halo. Our results cast doubt on the validity of predictions 
that long-lived resonances produce complicated structure in the Bose-Einstein correlation 
function at CERN SPS, when observed with a two-particle resolution of Qmin ~ 10 MeV. 

In the Appendix, we present one of the simplest possible application of the core-halo 
model to illustrate that the Gaussian "model independent" HBT radii are unreliable qual- 
itatively and are unreliable quantitatively in their frequently quoted original form, if the 
intercept parameter A < 1. At the end of the Appendix, we include the corrected definition 
of these parameters, that can be applied to core/halo type of systems with A < 1 also. 
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Appendix A: What do "Model-Independent Radii" measure? 

In this Appendix, we discuss the Gaussian model-independent radii of the HBT correla- 
tion functions in terms of the core/halo model, along the lines of ref. 

Let us recall the definition of the "model-independent" Gaussian radius parameters: 

Rh = {ix,-/34)ix,-/3,t))-{ix^-m){ixj-/3jt)), (22) 
C(K,Ak) = l + exp{-RljQiQj) , (23) 
iV(p) = (1), (24) 

{f{x,p)) = / d'^xf{x,p)S{x,p), (25) 



where i = side, out or long, (3i = Ki/Ko is the component of the four- velocity vector of the 
pair in the direction i and Qi = Aki, and S{x,p) is the emission function that characterizes 
a chaotic system, see ref. [^, [ll|] . 

For clarity, it is better to consider a one-dimensional system only. Let us consider a 
core/halo type of system, with Gaussian ansatz for both the core and the halo part, with a 
core radius of = 4 fm, a halo radius of Rh = 4:0 fm, and a core-fraction is assumed to 
be fc = 0.75. Let us compare the full correlation function C2{Q), given by eq. (|T^) to the 
core/halo model approximation C2^\Q), given as 

Cl'\Q) = l + X,eM-RlQ') (26) 
A, = fl = 0.56 (27) 
= = Aim (28) 

and to the model-independent Gaussian approximation to correlation function, C^((5), which 
reads as 

C^iQ) = l + Xcexpi-RlQ'), (29) 
Ac = 1 by definition, (30) 
i?^ = - ^ Rg = 20.3 fm. (31) 



This comparison is visualized in Figure |12. 



In Figure O it is clearly visible, that the model-independent Gaussian parameterization 
of the correlation function works only in the Q < h/Rh = 5 MeV region. The variance 
of the source function is dominated by the variance of the halo part, and this is picked up 
by the model-independent Gaussian radii. The approximation works in the unobservable 
Q range, however, in the well measurable Q > 10 MeV region the resulting correlation 
function approximates rather incorrectly the full correlation function. On the other hand, 
the core/halo model approximation fails in the unobservable Q < 5 MeV region, but it 
approximates with high accuracy the correlation function in the observable Q > 10 MeV. If 
a measurement is performed in the Q > 10 MeV range, the extrapolated correlation function 
will be a Gaussian with A < 1. We may conclude, that for a core/halo type of system, 
the Gaussian model-independent approximation yields qualitatively incorrect results in the 
range Q > %/ {x)\, when naively compared to a measurement of such a correlation function. 
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since A < 1 will be measured experimentally. If the radius parameter is considered only, 
its value is quantitatively wrong {Rg = 20.3 fm instead of 4 fm) in the model-independent 
approximation p^, since Rq = {x^) — (s)^ is dominated by the contribution of the halo. 



Note that the authors of refs. |jTT], ^ warned the readers against the straightforward 



applications of their result. Let us quote their warning from ref. [0]: " We would like to 
point out, that one must take care when comparing the above radii with the experimentally 
measured correlation radii, since the former measure second derivatives around Q = {], while 
the latter are parameters of a Gaussian fit to the whole correlation function ..." . 

Such non-Gaussian structures appear e.g. in the case of stable distributions, like the 
Lorentzian distribution, for which the Fourier-transformed distribution exists, however the 
first or the second moment of the distribution diverges thus the Fourier-transformed distri- 



bution is not analytic at small values of the relative momenta. See appendix of Ref. |jT6 
for further details on this point. Such a non-Gaussian correlation function may characterize 
even the core part of the distribution. 

From the above example it is clear that a direct application of the of the original defi- 
nitions of the Gaussian model-independent radii, as was proposed originally in refs. PB| , [TT 



leads to qualitatively and quantitatively unreliable results if applied even to very simple mod- 



els of core-halo type. The Gaussian model-independent radii were obtained in refs. p6| , |TT 
from an expansion around the Q = point. However, according to Figure |l|, in this region 
the correlation function may have a narrow and unresolved structure, dominated by the large 
regions of homogeneity in the halo. Obviously, the variances as determined by the above 
relationships shall be dominated by the variances of the halo part of the system, and in effect 
one obtains a nice Gaussian approximation to the Bose-Einstein correlation function — in 
the unresolved range of Q < Qmin ~ 10 MeV. However, the resolved part of the correlation 
function shall be missed completely by this approximation. See Figure for an illustration 
of the effect. 

From the above is also quite obvious, how to modify the original definitions of the model- 
independent radii to get an approximation to the measured correlation functions: the vari- 
ances should be evaluated for the core emission function only. Although the possibility of 
such a characterization was discussed in ref. ^ in a manner similar to the present Ap- 
pendix, the equations corresponding to such a description were not given there. We provide 
these below, since they may become a practically useful tool to characterize the two-particle 
Bose-Einstein correlation functions in high energy heavy ion collisions. 

This step results e.g. in a hydrodynamical like predictions for these radii and in a 
< A=K < 1 intercept parameter of the correlation function. A necessary condition for the 
applicability of the Gaussian approximation is that the first two moments of the source 
distributions of the core should be finite. The exceptional stable distributions (e.g. the 
Lorentzian) do not satisfy this criteria. We know that the mean of the core distribution is 
always finite hence the criteria is reduced to the requirement that the variances of the core be 
finite. Under these conditions, the two-particle correlation function in the core/halo picture 
can be rewritten as 

C'/\K, Ak) = 1 + A,(K) exp (-i?^ .(k) Ak.Ak,) , (32) 
A,(K) = [iVe(K)/Ar(K)]2, (33) 
RljiK) = {[Xi-l3it){x,-l3,t)),-{{xi-l3it)),{{x,-l3,t)),, (34) 
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(/(a;,p))c = J d^xf{x,p)S^{x,p), (35) 

where i = side, out or long as before, and Sc{x, p) is the emission function that characterizes 
the central core. Thus the halo contributes in this model to the reduction of the intercept 
parameter only and the variances of the core correspond to the Gaussian core/halo model 
radii of the measured correlation function. This result cannot be obtained with an expansion 
around Q = 0, on the other hand, it can be obtained with the help of a moment expansion of 
the source distribution function around x = {x)c, which becomes possible if the second mo- 
ment of the core distribution is finite. Hence, this result corresponds to a large Q expansion 
of the Bose-Einstein correlation function [^ . 

Although the above parameterization is a rather straight-forward combination of core - 
halo model and the expressions for the "model-independent" HBT radii, the results can still 
be considered simultaneously as a particular Gaussian approximation of the more general 
core/halo model result of eq. (n) as well as a generalization of the "model-independent" 
Gaussion HBT radii for systems of core/halo type. 
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Figure Captions 

Fig. 1 Full line indicates the effective, measured correlation function of the core/halo model 
using Gaussian ansatz for both the halo and the core, Rc = 4.0 fm/c, Rh = 40 fm/c. 
Dashed hne stands for the full correlation function, which includes the effect from the 
halo, resulting in a narrow and unresolvable peak if Qmin = 10 MeV is the experimental 
resolution. The extrapolated intercept parameter. A* thus deviates from the exact 
intercept of A^; = 1, however. A* carries important information about the fraction of 
core particles at a mean momentum. 

Fig. 2 The value of the three terms in the core-halo model correlation function, = 
0.75, Rc = 4.0 fm, Rh = 40. fm 

Fig. 3 The effect of on C2{Q) in the core-halo model for R^ = 4.0 fm and R^ = 40. fm. 

Fig. 4 The simulated A{Q) and B{Q) distributions, and the resulting C2{Q). Error bars are 
too small to be seen for the actual and the background distributions A{Q) and B{Q). 
Input values were Rc — 4.0 fm, Rh — 40. fm and fc — 0.75. To simulate the A{Q) 
distribution of the pion pairs, 349,476 points were sampled, for the B{Q) distribution, 
300,000. Note, that these A{Q) and B{Q) distributions were chosen to peak at 25 MeV, 
similarly to the CERN experiment NA44. 

Fig. 5 The behaviour of the fitted A* and R parameters as low Q data points are removed 
from C2{Q)- Solid line stands for A* = fc and Rfu = Rcore- 

Fig. 6 Radius parameters as a function of Qmin ^oi NA44's low pt data. 

Fig. 7 A* parameter as a function of Qmin for NA44's low pt vr+ data. 

Fig. 8 Radius parameters as a function of Qmin for NA44's high pt 7r+ data. 

Fig. 9 A* parameter as a function of Qmin for NA44's high pt tt"*" data. 

Fig. 10 Radius parameters as a function of Qmin for NA44's data. 

Fig. 11 A* parameter as a function of Qmin for NA44's data. 

Fig. 12 Comparision of the full correlation function (full line) to the core/halo model ap- 
proximation (dashed line) and to the model-independent Gaussian approximation (dot- 
ted line). The model- independent radii yield qualitatively and quantitatively unreliable 
results if A* < 1. 
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The influence of fc on the BECF 
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NA44 S+Pb - Low p.^ %* radius results 
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NA44 S+Pb - Low p-p intercept results 
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NA44 S+Pb - High p-^ %* radius results 
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NA44 S+Pb - High p.^ %* intercept results 
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NA44 S+Pb - radius results 
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NA44 S+Pb - intercept results 
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